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Abstract 

Q_i' We present a description of calculations of the electroweak amplitude for e + e — > tt 

process. The calculations are done within the OMS (on-mass-shell) renormalization 
scheme in two gauges: in R^, which allows an explicit control of gauge invar iance by 
examining cancellation of gauge parameters and search for gauge-invariant subsets of 
diagrams, and in the unitary gauge as a cross-check. The formulae we derived are 
realized in a FORTRAN code eef f Lib, which is being created within the framework of the 
project CalcPHEP. We present a comprehensive comparison between eef fLib results for 
the light top with corresponding results of the well-known program ZFITTER for the uu 
channel, as well as a preliminary comparison with results existing in the world literature. 
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Introduction 



The process e + e~ — > has already been studied for about ten years in connection with 
experiments at future linear colliders (see, for instance, the review []J). 

Actually, it is a six-fermion process (see 0); one of the channels is shown in Fig. [I]. 




Figure 1: The six-fermion e + e — > tt process. 

However, the cross-section of a hard subprocess, a(e + e~ — > tt), with tops on the mass shell 
is an ingredient in various approaches, such as DPA || or the so-called Modified Perturbation 
Theory (MPT), see §. 

In this article, we present a brief description of a calculation of the electroweak part of the 
amplitude of the e + e~ — > tt process. This calculation is a part of the project CalcPHEP H 
which started in Y2K after completion of the well-known project ZFITTER ||. One of main 
goals of this paper is to cross-check the CalcPHEP results against the results obtained with the 
other existing codes. 

As before, we use the OMS renormalization scheme, a complete presentation of which was 
recently made in J7J. However, for the first time we performed calculations in two gauges: 
and the unitary gauge. 

Note that there was wide experience of calculations in the gauge for processes such as 
H -> //, WW, ZZ, 7Z, 77, or e+e" -> ZH, WW. So, in @ and || a complete set of one-loop 
counterterms for the SM is given. Electromagnetic form factors for arbitrary £ are discussed 
1(J and [pjj . Explicit expressions can be found in the CERN library program EEWW [12 |. 
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However, we are not aware of the existence of calculations in the R^ gauge for the e + e — > tt 
process, although there are many studies in the £ = 1 gauge, see |jTB[ - [ |T6|j . 

Additional purposes of this study are: 

— to explicitly control gauge invariance in R^ by examining cancellation of gauge parame- 
ters, and search for gauge-invariant subsets of diagrams; 

— to offer a possibility to compare the results with those in the unitary gauge, as a cross- 
check; 
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— to present a self-contained list of results for one-loop amplitude in terms of Passarino- 
Veltman functions Aq, Bq, Cq and Do and their combinations in the spirit of the book |]7|, 
where the process e + e _ — > tt was not covered; this article may thus be considered as an 
Annex to this book; 

— to create a FORTRAN code for the calculation of the improved Born approximation (IBA) 
amplitudes and of the electroweak (EW) part of the cross-section of this process for a 
complementary study within the MPT framework; 

This article consists of five sections. 

In Section 1, we present the Born amplitude of the process, basically to introduce our 
notation and then define the basis in which the one-loop amplitude was calculated. We explain 
the splitting between QED and EW corrections and between 'dressed' 7 and Z exchanges. 

Section 2 contains explicit expressions for all the building blocks: self-energies, vertices 
and EW boxes. Note that no diagram was computed by hand. They are supplied by a new 
system, CalcPHEP, which is being created at the site brg.jinr.ru. It roots back to dozens of 
supporting form codes written by authors of [|?| while working on it. Later on, the idea came 
up to collect, order, unify and upgrade these codes up to the level of a 'computer system'. Its 
first phase will be described elsewhere [J5J. 

In Section 3, we describe the procedure of construction of the scalar form factors of the 
one-loop amplitudes out of the building blocks. One of the aims of this section is to create a 
frame for a subsequent realization of this procedure within the CalcPHEP project. 

Section 4 contains explicit expressions for the IBA cross-section. 

Finally, in Section 5 we present results of a comprehensive numerical comparison between 
eeffLib and ZFITTER. We also discuss some preliminary results of a comparison between 
eef f Lib and the other available codes. 
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1 Amplitudes 



1.1 Born amplitudes 



We begin with the Born amplitudes for the process e + {p + )e (p_) — > t(q_)t(q + ), which is 
described by the two Feynman diagrams with 7 and Z exchange. The Born amplitudes are: 



4? 



A\ 



eQ e eQtln ® 7m gi = -^ 47r «(°)^^7M ® 7m , 



1 



7„ 



Jf ) T+ - 2Q e s 2 w ® 7^ 4 3) 7+ - 2Q t s 



Q 2 + M| 



-ze 



if > J t (3) 7 M 7+ ® 7m7+ + <Uf V ® 7m7+ 



+4 3) 5 t 7 At 7+ ® 7^ + 4^7m ® 7m 
where j± — 1 ± 75 and the symbol ® is used in the following short-hand notation: 



;i.2) 



7m ( L i7+ + Qi) ® 7*> (^27+ + Q2) = v (p+) 7m ( L i7+ + Qi) « (p~) « (?-) 7*> (I27+ + Q2) « (?+) 5 

(1.3) 

furthermore 



-2Q/s 



2 



f = e,t 



;i-4) 



Introducing the II, QL, LQ, and structures, correspondingly (see last Eq. ( |1.2| )), we 
have five structures to which the complete Born amplitude may be reduced: one for the 7 
exchange amplitude and four for the Z exchange amplitude. 



1.2 One-loop amplitude for e + e 



it 



For the e + e~ — > ti process at one loop, it is possible to consider a gauge-invariant subset of 
electromagnetic corrections separately: QED vertices, 77 and Z7 boxes. Together with QED 
bremsstrahlung diagrams, it is free of infrared divergences. The contribution of QED diagrams 
is considered elsewhere [T7J. Here we keep in mind only the remaining one-loop diagrams 
forming electroweak corrections. The total electroweak amplitude is a sum of 'dressed' 7 and 
Z exchange amplitudes, plus the contribution from the weak box diagrams (WW and ZZ 
boxes) . 

Contrary to the Born amplitude, the one-loop amplitude may be parametrized by 6 form 
factors, a number equal to the number of independent helicity amplitudes for this process. 

We work in the so-called LQD basis, which naturally arises if the final-state fermion masses 
are not ignored []. Then the amplitude may be schematically represented as: 



i^+F* (s) + il»F e Q (s)} ® [z'7m7+^ (*) + *7m^q (*) + mJD^ (s) 



;i.5) 



If the initial-state masses were not ignored too, we would have ten independent helicity amplitudes, ten 
structures and ten scalar form factors. 
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with 

D li =(q+-q.) li . (1.6) 
Every form factor in the gauge could be represented as a sum of two terms: 

FU» 00 = F£l,» 00 + F% iD (s) . (1.7) 

The first term corresponds to the £ = 1 gauge and the second contains all £ dependences and 
vanishes for £ = 1 by construction. 

The LQD basis was found to be particularly convenient to explicitly demonstrate the 
cancellation of all ^-dependent terms. We checked the cancellation of these terms in several 
groups of diagrams separately: the so-called 7, Z, and H clusters, defined below; the W 
cluster together with the self-energies and the WW box; and the ZZ boxes. Therefore, for 
our process we found seven separately gauge-invariant subgroups of diagrams: three in the 
QED sector, and four in the EW sector. 

The 'dressed' 7 exchange amplitude is 

A™ A = i^f±a(s) lfl ®j,, (1.8) 

which is identical to the Born amplitude of Eq. ( |1.2|) modulo the replacement of a(0) by the 
running electromagnetic coupling a(s): 

OL 

a(s) = 7^- : ~ ■ (1.9) 



In the LQD basis the Z exchange amplitude has the following Born-like structure in terms 
of six (LL, QL, LQ, QQ, LD and QD) form factors: 

DC (s) [ 

A mA = • e 2 ——^ I 7/j7+ ® Jfl J+F LL (s,t)+ lfl ® 1»1+F QL ( S , t) 

+7 M 7+ ® 1^F LQ (s, t) + 7jU <g> 7 M F QQ (s, t) 

+7 M 7+ ® (-imtDr) F LD (s, t) + 7 M ® (-im*^) F QD (s, t) |, (1.10) 
where we introduce the notation for F { - (s, t) : 



F LL (s,t) = I^4 3) F LL (s,t), 

F QL (s,t) = 5 e 4 3) F QL (s,t), 

F LQ (s,t) = iPS t F LQ (s,t), 

F Q Q ( S ^) = $e5 t F QQ (s,t), 

F LD (s,t) = I®4 3) F LD (s,t), 

F QD (s,t) = 6 e 4 3) F QD (s,t). (1.11) 
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Note that tilded form factors absorb couplings, which leads to a compactification of formulae 
for the amplitude and IBA cross-section, while explicit expressions will be given for untilded 
quantities. The representation of Eq. ( |1.10|) is very convenient for the subsequent discussion 
of one-loop amplitudes. 

Furthermore, in Eq. ( l.lOj ) we use the Z/7 propagator ratio with an s-dependent (or 
constant) Z width: 



Xz{s) 



4q2 r 2 



;i.i2) 
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2 Building Blocks in the OMS Approach 



We start our discussion by presenting various building blocks, used to construct the one-loop 
form factors of the processes e + e~ — > ff in terms of the A , B , Co and Do functions. They 
are shown in order of increasing complexity: self-energies, vertices, and boxes. 

2.1 Bosonic self-energies 

2.1.1 Z, 7 bosonic self-energies and Z— 7 transition 

In the i?£ gauge there are 14 diagrams that contribute to the total Z and 7 bosonic self-energies 
and to the Z—j transition. They are shown in Fig. 0. 

With S zz , S Z1 and S" 77 standing for the sum of all diagrams, depicted by a grey circle in 
Fig. 13, we define the three corresponding self-energy functions T, AB : 



Szz — 


(277 


\/ 2 2 S 
lOTT^C^ 


Szj 


(277 


167T 2 C W 


7 = 


(277 


1 c 2 77 



ZZ ) 



Z7 j 



(2.1) 
(2.2) 
(2.3) 

All bosonic self-energies and transitions may be naturally split into bosonic and fermionic 
components. 

• Bosonic components of Z, 7 self-energies and Z-j transitions (see diagrams Fig. 0) 



Z \3R Z V2 



c z -9c 4 



(2.4) 



(1 



2c 4 . 



6 



J zz 1 



.) = ^ 

; 12 
+ 



Ac 2 w (5 - 84 - 12ct) + (l - 4c 2 , - 364) TT 



— + 10 - 2r HZ + (r HZ - l) 2 i2 z 



B*{-8;M B ,M Z ) 



B^(-s;M w ,M w ) 



+ 



+ 1 + (1 - r HZ ) R z 



r. 



7-(l-r r 



L,(M 2 H ) 



+2c 



l + 8c 2 w -24c 4 w )L ll (M WJ 



(2.5) 



+ -(1 -2c 2 ,) — -6(l + 2ct 



- 3(1 + 2c 2 ,) - 9r flZ - (1 - r HZ ) 2 
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Z,7 Z,i 




u,d W + Z 

\NWb ^A/VW + WW* WW/W + WW# WWW 

(1) u,3 (2) W~ (3) if 

0+ 

+ WW# www + ww» www 

(4) 0" (5) W" 

0° 0+ 

+ vwW www + vw» www 

(6) If (7) 1 

+ WW«i Www + VWVli i»AA/W 

(8) X- (9) X + 



(10) W (11) H 




+ VWWVWWVWW + WWWMVWWW 



(12) 0+ (13) 

/ \ 

I I 



o 



+ VWWVA^WWA/W + VWWVMVWWW 



( 14 ) A (Z) 
+ VWWWVWA/WW 



Figure 2: (Z, 7)-boson self-energy; Z-7 transition. 
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Here L M (M 2 ) denotes the log containing the 't Hooft scale \i: 

M 2 



In 



(2.6) 



and it should be understood that, contrary to the one used in [0, we define here 

B ( -s;M ll M 2 ) = - + <( - s; M x , M 2 ), 



(2.7) 



meaning that Bq also depends on the scale \i. We will not explicitly maintain /i in the 
arguments list of and Bq. Leaving /z unfixed, we retain an opportunity to control /i 
independence (and therefore UV finiteness) in numerical realization of one-loop form factors, 
providing thereby an additional cross-check. 

Next, it is convenient to introduce the dimensionless quantities II z ° s (,s) and Il^° s (s) (vac- 
uum polarizations): 



S b z ° s ( S ) 



-*n b / 7 s ( s ), 



= sll^(s). 

In Eqs. ( ^.5| ) and ( ^.7[ ) and below, the following abbreviations are used: 



M 2 
Ml ' 



r ..-< R -Ml 

rrij s 



Rp. 



Ml 



(2.8) 
(2.9) 

(2.10) 



Since only finite parts will contribute to resulting expressions for physical amplitudes, which 
should be free from ultraviolet poles, it is convenient to split every divergent function into 
singular and finite parts: 



n b ° s - F ( s ) 



3- + n bos < F ' 



77 



s), 



(3 + AR W ) B%( - s - M w , M w ) + AR W L^M, 



and 



n b ° s (s) 



Z7 



l - + 3c 2 w + 2ic) \ + n h ™' F (s) 

I + 3cL + 4 (I + cl) R w 



1 
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Bi 



s;M w ,M w ) 



Ac 2 w )R w L^M w) 



With the Z boson self-energy, E zz , we construct a useful ratio: 

1 S zz (s) - S zz (M 2 Z ) 



Ml - s 



(2.11) 
(2.12) 

(2.13) 
(2.14) 

(2.15) 
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which also has bosonic and fermionic parts. The bosonic component is: 



V h ™ (a) 



V 



bos.F 



— (-- : + -ci +3r,' ( ]- + P 

1 



^bos,F 



— + -c z 
12 3 w 



— c 4 -4c 6 



(2.16) 
(2.17) 



Ml 



x 



M| - s 



B£( - s; M w , M w ) - B£( - M 2 Z - M w , M M 



1 



rC,, 



+ 1 



12 3 
1 



3c 4 w )bS(-s;M w ,M w ) 



-J HZ + 



1 

12 



[l-r 



— r 

12 HZ 

2 



Ml 



M 2 Z - s 



B${ - s; M H , M z ) - B£( - M\- M H , M z ) 



Bn 



s; M H , M z 



~R M (1 - r HZ ) [r HZ (L,(Ml) - l) - L,(M 2 Z ) + 1 
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1 - 2d 



• Fermionic components of the Z and 7 bosonic self-energies and of the Z-7 transition 

These are represented as sums over all fermions of the theory, They, of course, depend on 
vector and axial couplings of fermions to the Z boson, Vf and a/, and to the photon, electric 
charge eQf, as well as on the colour factor c/ and fermion mass m/. The couplings are defined 
as usual: 



v f = if - 2Q f s 



2 

W 5 



a f 



r(3) 
'/ ' 



with weak isospin if, and 



2 1 
Qf = —1 for leptons, +- for up quarks, — - for down quarks, 

o o 

Cf = 1 for leptons, 3 for quarks. 
The three main self-energy functions are: 



/ 



(v 2 + aj^j sBf( — s;mf,mf) — 2a 2 j : m'jBo( — s; mf, mf) 



(2.18) 

(2.19) 
(2.20) 

(2.21) 

(2.22) 
(2.23) 



The quantities IT^ and IFj* are different according to different couplings, but proportional 



to one function Bf (see Eq. (5.252) of for its definition): 

= AY^c f Q)B f {-s ]m} ,m f ) 



77 

n fer ( 



2 J2 c fQf v f B f (-s;m f , 



,m f ) 



[2.24) 
[2.25) 
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As usual, we subdivided them into singular and finite parts: 
n» = -l(^N f -4slJ2c f Q^-_+U%' F (s), 



(2.26) 



--s 2 w )N f + 4siJ2 c fQ 



In Eq. ( |2.26|) , Nf = 24 is the total number of fermions in the SM. We do not show explicit 
expressions for finite parts, marked with superscript F, because these might be trivially derived 
from Eq. ( 2.21 ) and Eqs. ( 2.24j ), ( 2.25| ) by replacing complete expressions for Bj and B with 
their finite parts B F and Bq , correspondingly. 

2.1.2 W boson self-energy 

Next we consider the W boson self-energy, which is described, in the gauge, by 16 diagrams, 
shown in Fig. |3|. 

First, we present an explicit expression for its bosonic component: 



S bos (s) = M. 

WW V ' 1 

where 

E bo S ,F ( 



19 1 1 



6 R w 4 



Ml 



6 



+ 2 



- + 10 + 3r HW 



+ S bos ' F (s) 

WW v J 



12 



l-40d) — + 2 — -27-8c 2 v 

■tiw \ C w 



(2.27) 



+ - 



+ 8\R V 



B ( -s;M w ,M : 



+ 



R\ 



+ 2(h-r HW ) + (l-r HW fR w 



B (-s;M w ,M B 



-8si 



5 

Rw 



+ 



HW 
1 

2~ 



+ 2-RS) B (-s;M w ,0) 
L,{M 2 H ) 



7- 0--r HW )R w 



18 1 , , / 1 

— r + 1 - 164, + 4 — 



8 



+ 



4 1 

3 Rw 

1 



3 — + 2 - 9r f 



+ 6s 



1 -n 



^HW ''fiff; 



(2.28) 



12 



u 



W+ 



W+ 



W 1 VV / -"~\ ^% 

VWW jWWV = VWV* WVWV + VWV* «AMA/ + VW\A# MAM/V 

A* 7 f ^Al^ 

(1) 3 (2) Z (3) 7 



</> H 



<t>" 



+ VWVw WVWV + VWVfc, AAAM/ + VW\AA AAMA/ 



(4) If 



(5) Z 



(6) 7 



(7) # (8) 0° 



^,7 ^ + 



+ VW\Af *A/WV + WW* iWWV, 

(9) X- (10) C 



'11) W 



12) Z 



+ 




+ 




(13) 7 



+ 




(14) if (15) 0+ (16) 0° 

li ii il 

+ VW\AAA^WWWV + VWVW^WWm + V¥VWV\^WWWV 



(17) A 

+ ^A/W\AA%W\AAAA 



Figure 3: boson self-energy. 
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Secondly, we give its fermionic component: 



-i fcr 



(s) = s^CfBfi- s;m f/ ,m f ) + ^c / mJ J B 1 (- s;m r ,m f ), (2.29) 

f=d f 



where summation in the first term extends to all doublets of the SM. 
2.1.3 Bosonic self-energies and counterterms 

Bosonic self-energies and transitions enter one-loop amplitudes either directly through the 
functions V z (s), U. Js) and H z (s), or by means of bosonic counterterms, which are made 
of self-energy functions at zero argument, owing to electric charge renormalization, or at 
p 2 = —M 2 , that is on a mass shell, owing to on-mass-shell renormalization (OMS scheme). 

• Electric charge renormalization 

The electric charge renormalization introduces the quantity z 1 — 1: 

2 



> 7 - 1; 



with bosonic (see Eq. (6.161) of 0): 

/ -| \bos 

[z 1 - 1) 

and fermionic 



n 77 (o) 



M 2 



S 3Q (0) 



(2.30) 



3[-_-L,(M 2 w : 



, fcr 



(2.31) 



^2.32) 



components. 

• p-parameter 

Finally, two self-energy functions enter Veltman's parameter Ap, a gauge-invariant combi- 
nation of self-energies, which naturally appears in the one-loop calculations: 



Ap 



Ml 



,(M 2 W )-E ZZ (M 2 Z ) 



(2.33) 



with individual components where we explicitly show the pole parts: 

/ 1 



bos 



Ap 

Ap fcr 



f + 74.)i + Ap b ^, 



fer.F 



(2.34) 

(2.35) 
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The finite part of Ap bos is given explicitly by 
/ 1 4 17 



Ap 



bos,F 



12c w 3c^ 



B* ( - Af£; M w , M z ) - c 2 w B* ( - M 2 ; M w , M w ) 



+ ( 1 " \r HW + irj^) <( - M w , M ff ) 
-II- gr« + y^L) 7£ 



( - M z \ M z , M H ) - 4s 2 w B* ( - M w , 0) 



+ 



12 



- 24 + r HW L,{Ml) + L M (M^) - 1 



( 1 9 , \ , 9N 1 19 22 

- — + 14 + 1&& - 484 + r HW L,(M 2 W ) - — -— + — 



C- 
1 i 



3c- 



^2.36) 



while the finite part of Ap fcr is not shown, since it is trivially derived from the defining equation 
fl2.33|) by replacing the total self-energies with their finite parts. 

2.2 Fermionic self-energies 
2.2.1 Fermionic self-energy diagrams 

The total self-energy function of a fermion in the gauge is described by the six diagrams 
of Fig. |. 



X 




X 



X 




X 



1) 7 



+ 



-X/\x- 



(4) H 




+ 



(2) Z 

^x\^ 

(5) 0° 



X 




X 



(3) W 

+ ^rvx 

(6) 



Figure 4: Fermionic self-energy diagrams. 



Calculating derivatives straightforwardly and substituting the a^'s, we obtain (see 
explicit expressions for the wave-function renormalization factor yfz7~^. 
It is convenient to distinguish the electromagnetic components 

,2 \ 



(\ < in i , < ni ., ., ; 1 13 vnfi 



and the weak components 



- I ={w v ± w a ) , 



[2.37) 



[2.38) 
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where 



u>. 



w, 



w„ 



1 1 



.__l(„ ?+a?+2a?rlz )I 



£ 



+2 (1 + 2r tz ) Mf 



m\\m u M z ) - L^m 2 ) 



— (<( - m 2 t ; M z , m t ) + L M (M 2 ) - 
,(m 2 )] + 2a 2 t r tz \— (<( - m 2 ; M z , m t ) 



1 

~16 
+2(1 



t (M 2 ) - l) - 6M 2 z B 0p ( - m 2 ; m u M z ) - L ' 2 
(2 + ^)U + 1+r - 



e 



r tw 

-r tw )MlB 0p (-m 2 t ;0,M w ) + 



'B£{-m 2 ;M w ,0) 



+ 2a 2 t 

'Ml) + 1 
i 



Vrtz v 



e 



1 1 

~T^ r ^\ - + 
16 I £ 

-2 (4r tj 



? F ( - m 2 ; M H , mt) + L^(M 2 H ) - I 



'tw 

1) M 2 H B 0p ( - m 2 ; m t , M H ) - L^m 2 ) + 1 
~u t a t { - - — [<( - m 2 ; M z ,m t ) + L^(M 2 ) • 



+2B F (-m 2 ;M ; 



z ,m t ) + L^m 2 ) - 2 



3^ [B^(-m 2 ;M w ,0) - 



1 



-r tw B^(-m 2 ; M w ,0) 



\r tw 



1 ) ^« 



(2.39) 



(2.40) 



(2.41) 



(2.42) 



(2.43) 



2.3 The Zff and 7// vertices 

Consider now the sum of all vertices and corresponding counterterms whose contribution 
originates from the fermionic self-energy diagrams of Fig. [|. This sum is shown in Fig. |5|. 





+ 




Figure 5: Zff and 7// vertices with fermionic counterterms. 
The formulae which determine the counterterms are: 

2(v^-/), 



Q 



(2.44) 
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FT' = (V^-J)-(V^-/), (2.45) 

= 5}(^-I), (2.46) 

F z ^ = a 2 f (^-I)-5 2 f (^-I), (2.47) 

where 

6f = v f - a f , a f =Vf + a f . (2.48) 

For the sum of all 7 — > // and Z — > // vertices (the total j(Z)ff vertex depicted by a grey 
circle in Fig. |5|) we use the standard normalization 



m 2 = (27r) 4 ^, (2.49) 



and define 



V](s) = (2*)S'— G M ( S ), (2.50) 
V*(s) = (2n)U-^-Z(s), (2.51) 



while we denote the individual vertices as follows: 

G>(a) = (s) + (s) + Gjf (s) + G^ (s) , (2.52) 

V s ) = Z;(s) + Z*(s) + ZW(s) + Z*(s). (2.53) 

All vertices have three components in our LQD basis. 

2.3.1 Scalar form factors 

Now we construct the 24 = (4 : A, Z, H, W -virtual) ® (3 : L, Q, D) ® (2 : 7, Z -incoming) scalar 
form factors, originating from the diagrams of Fig. [|. They are derived from the following six 
equations — three projections for 7// vertices: 

F2 B {s) = -^{G°(s)[ig 3 7 , 1+ ] + s w Q f F2> ct }, (2-54) 
s w l j 

F^is) = -^{Gl{s)[ig"l»\+s w Q f F%<*}, (2.55) 
FT(s) = -^GZWtfmtlDr], (2.56) 
and three projections for Zff vertices: 

Fl B {s) = ^\z*(s)^ 1 , 1+ } + -fA, (2.57) 
l f { c w ) 



F z Q B {s) = ^f{ ZB A^[i9 3 l,) + ^F^l (2.58) 
F?{8) = ^Z»(s)[g 3 m t ID,}. (2.59) 
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Here we have / = t,e, and B = A, Z, W, H, and we introduce the symbol [. . .] for the definition 
of the procedure of the projection of G (s) and (s) to our basis. It has the same meaning as 
in form language [|(J, namely, e.g. (s) [ig 3 ^^ + ] means that only the coefficient of [ig 3r y^ + ] 



of the whole expression G^ (s) is taken (projected). 

The factors l/(QfS w ), 2/{s w l[ ) and 2c^ v /(s w I^) for 7// vertices, and the factors 
2c w /It 3 \ 2c w /5f and 2c w jlf ) for Zff are due to the form factor definitions of Eq. ( |1.10| ). 



The total "ftt and Ztt form factors are sums over three bosonic contributions B = Z,W, H 
since we separated out the contribution of the diagram with virtual 7 = A: 

F2 tt {s) = FT (s) + F2 W {s) , 

F^{s) = F^ z (s) + F^ w (s) + F^ H (s) , 

F2 a {s) = F2 Z (s) + FT (s) + F2 H (s) , 

Ff{s) = F? (s) + FT (s) + Ff (s) , 

Ff ( S ) = F?(8) + F™(8)+F?(a), 

F2 u (s) = F2 Z (s) + F2 W (s) + F2 H (s) . (2.60) 
The quantities F2&,d (s) originate from groups of diagrams, which we will call clusters. 



2.4 Library of form factors for Btt clusters 

Here we present a complete collection of scalar form factors F2^2> (s) originating from a vertex 
diagram with a virtual vector boson, contribution of a scalar partner of this vector boson, and 
relevant counterterms. 

Actually three gauge-invariant subsets of diagrams of this kind, A, Z and H, appear in 
our calculation. They may be termed clusters, since they are natural building blocks of the 
complete scalar form factors, which are the aim of our calculation. Again, in the spirit of our 
presentation, we write down their pole and finite parts. The remaining vertices with virtual 
W and <f) + , cf)~~ with relevant counterterms we also define as the W cluster. However, the latter 
diagrams do not form a gauge- invariant subset. 



2.4.1 Form factors of the Z cluster 

The diagrams shown in Fig. |^ contribute to the Z cluster. 

Separating out pole contributions 1/e, we define finite (calligraphic) quantities. We note 
that, if a form factor F l J(s) has a pole, then the corresponding finite part J-^{ s ) is /^-dependent: 

F2 z (s) = T2 z {s), 

F2 z (s) = n z {s), 

F^s) = 
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X 




X 



+ 



0° 



Figure 6: Z cluster. The two fermionic self-energy diagrams in the second row give rise to the 
counterterm contribution depicted by the solid cross in the last diagram of the first row. 



Here the finite parts are: 



1 



Qtv t 2 2 



1 



M 2 Z C Q ( - m 2 , -m 2 , -s; m t , M z , m t ) 



(2.61) 



(2.62) 



-3B F { - s- m t , mt) + 2<( - m 2 ; m t , M z ) - L^m 2 ) 

M 2 

+B^(-m 2 ;m t ,M z )-2(l + 4r tz ) ' 



z L ab (m t ,m t ,M : 



ZJ (1 



2 (2 (1 - r tz ) + M 2 Z C ( - m 2 , -m 2 , -s; m t , M z , m t ) 
B F { - s- m t , mt) + 45 F ( - m 2 ; m t , M z ) + L^{M 2 Z ) - r tz B F 2 { - m 2 ; m t , 
(1 + 2r tz ) M 2 z B 0p ( - m 2 ; m t , M z ) - ^ 



*3r 



) M 2 Z C Q ( - m 2 , -m 2 , -s; m t , M z , m t ) 
-3B F ( - s; m t , m t ) + 4B F ( - m 2 ; 



M z ) 



+2v t a t r tz -4M 2 Z C ( - m 2 , -1 



m 2 t1 -ml -s; m t , M z , m t ) + 2 (L M K) - L^M 2 Z )) 

■2(l-r tz )Bg(-m 2 ;mt,M z ) + l 
9 / 



-2 (1 - r tz ) B r d2 ( - mt; m u M z ) + 1 

— ( (1 + 2r tz ) M 2 z B 0p ( - m 2 ; m t , M z ) + \ 
r tz V * 

+2a 2 r tz B F ( - s; m u m t ) + L^M 2 ) - r tz B F d2 { - m 2 ; m u M z ) 
+6M 2 B 0p (-ml,m t ,M z )-^ 



-4 



1 M 2 1 

(Av t a t - afj r tz ^L ab (m t , m t , M z ) \, 



(2.63) 
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- AM 2 z Cq( - m 2 t , -ml, -s; m t , M z , m t ) 



2Q t 1 \ v 2 t +a 2 t 
I^cl^A 2 |_ 
+5 F ( - s; m t , m t ) ~ ^( - m 2 t ; m t , M z ) - L 
M 2 

+2 + 6-^L ab (m t ,m t ,M z ) 

+a 2 t 2 (3r tz - — ^ M 2 C ( - m\, -m 2 t , -s; m t , M z , m t ) 
+B?( - m 2 t ; m t , M z ) + L„(M|) - 1 - 

(Tn 2 \ ' 
2-3^-jL a6 {m t ,m t ,M z 



^( m t) + B di(-m 2 t ;m t ,M z ) 



Bq( - s; m t , m t ) + L^{m 2 t ) - 2 



(2.64) 



1 <3v 2 + 



a; 



4^1 



2 ^2 + i-) - 2r^ M 2 C ( - m t 2 , -m 2 , - s; m t , M z , m t ) 
s; m t , m t ) + 8Bq( - m 2 ; ra t , M z ) - L^(m 2 t ) 



— m 2 



m 2 t ;m t ,M z )-2 



+2B F 



2 ;m t ,M z )-2(l + 2r tz )M 2 z B 0p ( 
4m 2 C ( - ml, -ml, -s; m t , M z , m t ) + 3r 
^; m t , M z ) + 3L M (M 2 ) + 2L M (m 2 ) 
I; "it, M z ) + 2 (1 - 7r tz ) M 2 z B 0p ( - m 2 t ; m t , M z ) - 



Bq(- s;m t , m t ) + L^m 2 ) 



■ml; 



— m 2 



'di 

M 2 

2 [(?>v 2 t + a 2 ) (1 + 4r tz ) - 2a 2 t r tz ] —^L ab (m t , m t , M, 



z) , 



(2.65) 



T z Q z {s) = J-{j<S t 2 2(2- 2r tz + -^)M 2 z C (- m 2 t ,-m 2 t ,-s;m t ,M z , m t ) 

-3B£ ( - S ; mt, mt) + 4<( - m 2 t ; m t , M z ) + L^m 2 t ) - B£( - m 2 t ; m t , M z ) 

M 2 

-2 (1 + 2r tz ) M 2 z B 0p { - m 2 t ; m t , M z )-2- 2-^L ab (m t , m t , M z 

^3r 



+a t r tz vt 



1 

+ — 

r tz 



B' 



-2M 2 C ( - ml, -ml, -s; m t , M z , m t ) 



ml; m t , M z ) + L^m 1 ) -!-£&(- m 2 t ; m t , M z ) 



- (1 + 2r tz ) M 2 B 0p ( - ml; m t , M z ) 



M 2 

Q-r^L ab (m t ,m t ,M z ) 
A 3r 
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--at 
2 



8Mf C ( - m 2 , -m 2 , -s; m t , M z , m t ) 



-B%{ - s; m t , m t ) - L^m 2 ) + 2 + B£( - m 2 ; m t , M z ) 

M 2 

-QM 2 z B 0p ( - m 2 ; mt, M z ) - 10-^L afe (m t , m t) M z 

^3r 



5, 



4M 2 C ( - ra 2 , -^ 2 , -s; m t , M z , m t ) 



Ml 



-Bq( - s; m t , m t ) + 1 - 6^--^L afe (m t , m i; M z 

^3r 



(2.66) 



— m 7— I ( 3 fl2 + v t 



-4M 2 C ( - m 2 , -m 2 , -s; m t , M z , m t ) 



t 

+5 F ( - s; m t , mt) ~ 25 F ( - m 2 ; m*, M z ) - L M (m 2 ) + - m 2 ; m t , M z 
M 2 

+2 + Q—^L ab (m t , m t , M z 

^3r 



+2v t a 2 



2 (Vr iz - 2-^-) M 2 C ( - m 2 , -m 2 , -s; m u M z , m t ) 



+B%( - m 2 ; m t , M z ) + L^(M 2 ) - 1 - r u 



-2(4-3^jL afe {m t ,m t ,M z 



Bjf ( - s; m t , m t ) + L^(m 2 ) - 2 



(2.67) 



In Eq. ( |2.63| ) the 'once and twice subtracted' functions B^ and B^ 2 are met: 



B Ff ™2 



'til* 



m t ;m t ,M z ) = — 



B^{-m 2 ;m t ,M z ) 



Ttz 
1 

r 2 tz 



<( - m 2 ; m t , M z ) + L M (M 2 ) - 1 
B*{ - m 2 ; m t , M z ) + L M (M 2 ) - 1 



-r tz [L,(m 2 )-L IM (M 2 ) + 



(2.68) 



1 

2 

They remain finite in the limit m t — > 0. 

We note that, for the Z cluster, all the six scalar form factors F2^d (s) are separately 
gauge-invariant . 
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2.4.2 Form factors of the H cluster 

The diagrams of Fig. [j] contribute to the H cluster, 




H + 




/ + 




{Z)H/ 




H 



Figure 7: if cluster: the vertices and the counterterm. 



Separating UV poles, we have: 



F Z Q H (s) 



11 11 



T (s) , 



lG\Qt 



r tw ~~ r~ ~ + ^ H l s ) 



with the finite parts: 



1 



Fq H ( s ) = o r tw{ 8m t C o(- m t>- m t,- s ', m t,M H ,m t ) 



+Bq( - s; m t , m t ) + L^(m 2 t ) - 2-B$ 1 (- m 2 ; m t , M H ) 

M 2 

-2 (1 - Ar tH ) M 2 H B 0p ( - m 2 ; m t , M H ) - 2-^L ab (m t , m u M E 

^3r 



22 



FZ H (s) = -%^{-6MlC (-m 2 t ,-m 2 t ,-s;m t ,M H , 

2r t ' ^3r { 

-3Rf ( - .<?: m*. - AR? ( - m 2 : m.+ . M„) - J, 



+3Bq( - s; m t , m t ) - 4B,f ( - m\; 

M 2 ) 

- m 2 ; m t , M H ) + 6-^L ab (m t , m t , M H ) , 

.F£ H (s) = K tw ^Am 2 C (-m 2 ,-m 2 ,-s;m t ,M H ,m t ) 

+ [4(1- r tz ) + (1 - r^) 2 ^] M 2 Z C ( - m 2 , -m 2 , -s; M H , rr 



m t ) 

2 ;m t ,M H )-L^m 2 ) + 2 



\{l-r tz ) + {l-r HZ ) 2 R 



M 2 C ( - m 2 , -m 2 , -s; M H , m t , M z ) 



+2<( - s; M z , M H ) - i<( - s; m t , m t ) + h,(m 2 ) 



+ {l-r HZ )R z 



M. 



(2.71) 



- m 2 ; m t , M H )-{1- Ar tH )M 2 H B Qp ( - m 2 ; m t , M H ) + 2 
'<( - m 2 ; M z , m t ) - <( - m 2 ; m t , M H 

{r HZ - 8r tz ) L ab (m t , m t , M H ) - 2 (3 - r HZ + 4r tz ) L Hi (m t , M H , M, 

^3r L 



(2.72) 



? Z q ( s ) = r tw {m 2 t C Q {- m 2 t ,-m 2 t ,-s;m t ,M H ,m t ) 



+ 



+- 



1 + (1 - r HZ ) R z ] M 2 Z C ( - m 2 , -m 2 , -s; M H , m t , M z ) 
Bq(-s; m t , m t ) + L^(m 2 ) - 2 



+R Z 



B%{ - m 2 ; M z , m t ) - <( - m 2 ; m t , M H )] - -££( - m 2 ; m t , M H ) 



1 1 M 2 

-7 (1 - 4r to ) M|fi 0p ( - m t 2 ; m t , M H ) - --r^L ab (m t , m t , M H ) 

4 4 A 3r 



+ 



la t 



4 5, 



4r tz + (3 + r ffZ ) (1 - r HZ )R z \ M 2 Z C ( - m 2 , -m 2 , -s; M H , m t , M z ) 

+<( - s; m t , mt) - 2<( - s; M z , M H ) - 3 
+ (3 + r HZ ) R z \b£( - m 2 ; M z , m t ) - B^( - m 2 ; m t , M H 



(2.73) 



-2- 



Ml 



*3r 



(m t , m t , M H ) - (3 - r HZ + 4r tz ) L Hi (m t , M H , M z 



2/f ) c 2 v A 3r 



-3r tz M|C ( - m 2 , -m 2 , -s; m t , M H , ra t ) 



+2 



2(r, 



1) — - r ffZ + 2r tz 
s 



M 2 C ( - m 2 , -m 2 , -s; M H , m t , M z ) 
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2 

—4 — - 

s 



£ F ( - m?; M z , mt) ~ - m 2 - m t , M H ) 



2B^-s;M z ,M H ) 



-2B^(-m 2 ;M z ,m t ) + -r tz 



Bq{ - s; m t , mt) ~ - m 2 ; m t , M H ) 



+ -r HZ 



B%{ - m 2 - m t , M H ) + L M (M 



B 



m 2 ;m t ,M H ) + L^(m 2 ] 



M 2 1 

+ 3r tz -^L ab (m t ,m t ,M H )\. (2.74) 



Again, the five (one does not exist) scalar form factors in Eq. ( |2.69|) are separately gauge- 
invariant. Note also that UV poles persisting in the scalar form factors of the H cluster cancel 
exactly the corresponding poles of the Z cluster. In other words, the form factors of the 
'neutral sector' cluster [Z + H) are UV finite. 

In total, we have 11 separately gauge- invariant building blocks that originate from Z and 
H clusters. 



2.4.3 Form factors of the W cluster 

Finally, the W cluster is made of the diagrams shown in Fig. R 




■r + 




+ 




/' + 




+ 




+ 



+ 



/' + ww* if + 




Figure 8: W cluster: the first row shows the abelian diagrams of the cluster, the last row the 
non-abelian diagrams; the second row shows diagrams that contribute to both counterterm 
crosses (last diagrams in first and third rows). 

In the formulae below, we present contributions to scalar form factors from all the diagrams 
of the W cluster, not subdividing them into abelian and non-abelian contributions. To some 
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extent two sub-clusters are automatically marked by the type of arguments of Co functions 
and typical coupling constants. Separating poles, we have: 



21, 



(3) 



/ x2 2 

3+(l + r 5 ) + — 



■o(6 + r«) 



- 2-r 



M W C ( - m 2 , -m 2 , -s; m 6 , M W) m b ) 
Bq( - s; m b , m b ) - Bq{ - m 2 ; m b , M w ) 
B dw( - m t ; M w) - 1 



+ 



4 - (2 + r"J (3 + 3r tw + r bv 



^r^L ab (m t , m b , M M 

^3r 



-r 6w (2 - r tb ) M W C ( - m 2 , -m 2 , -s; M w , m b , M w ) 
+\ (6 - r a 



Bg( - s; M w , M w ) - Bl{ - m 2 ; m b , M w 



+2B* ( - mi; m b , M w ) + - (2 - 3r£ 



B dw{ - m ti m b, M w ) + 1 



i 4 - 

Qt 



41, 



(3) 



where 



2 + r tw (13 + r tw ) - r bw (1 + r 6w ) 
Ajr Im - z\ilmB$( - m 2 ;M w ,m b ), 

Af Im = (1 - r bw ) h r tw . 



Ml 



A 



3r 



!L„ a (m t ,m b ,M„ 



(2.75) 



(2.76) 



The last term in Eq. ( p. 75 ) is due to a non-cancellation of the imaginary part of the function 
Bq( — m 2 ; M w ,m b ) which appear in real counterterms and complex- valued vertices. 



4Q 



-4M W C ( - m 2 , -m 2 , -s; m b , M w , m b ) 



t { 



+S^( - s; m&, m b ) - Bq( - m 2 ; m b , M w ) - 1 



n 



/ (1 - tw) (2 + r bw ) i _ 



M 



r 4w + 2r blv M WJ B 0p ( - m t 2 ; m 6 , M w ) 



+2 (3 + r~ b ) -r^Lah (m t , m 6) M M 
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(3) 



A 3r 

,2^y / 2 2 



2m b C ( - m t , — ra 4 , — s; M w , m b , M w 



-5 F ( - s; M w , M w ) + <( - m 2 ; m 6 , M w ) - 1 + 



?7 



+ 



'2 - r bw , (1 + r fev 
% 



1 - r tw + 2r bH , M WJ B 0p ( - m 2 ; m 6 , M M 
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+2 (7 + r+J (m t , m 6 , Af w ) 



+ ^A^i ImE F ( - m^; M w , m b ) , (2.77) 



i f Q 6 



A 3r 1 2L (3) 



-2 



_2 2 
8 + r tb - Qr bw + — 

ri w 



+ (10 + r tw - 3r bw ) 



Bq( ~ s; m b , m b ) - B£( - m 2 ; m b , M w ) 



+ (2 + rj) 



B^ w (-m 2 ;m b ,M w ) + l 



(l + r tw )(2 + r t 

w ) r bw (l + r b w J 



+6 
+2 

xMlC {-m 2 t ,-m 2 ,-s;M w ,m b ,M w ) 



M 2 

-r^L ab (m t ,m b ,M w ) 

A 3r , 



(l - r^) (2 - n w ) - r 6w - 4r 6w - 



+ (6 - r tw - 5rw) 



B£( - s; M w , M w ) - B^( - m 2 ; m b , M w ) 



+ (2 + rt b ) 



Bdw(- m t', m b,M w ) - 1 



(1 - r tw ) (2 + r tw ) ~ r bw (1 + 2r tw + r bw ) 



M 2 

X yJ - T* 



L na (m t ,m b ,M w )\, (2.78) 



*3r 
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(3) 



/ \ 2 2 

3+(l + r») 



"2 ( 6 + ^ 



M^C ( - m 2 , -m 2 , -s; m b , M w , m b ) 
Bq( - s; m b , m b ) - Bq( - m 2 ; m b , M w ) 



+ 2 ( 2 " r tb) [ B L( ~ m 2 t ; m b , M w ) - 1 
+ 4 - (2 + r- b ) (3 + 3r tw + r bw 



■L ab (m t ,m b ,M„ 



+r bw M' w C ( - m 2 , -m 2 , -s; m b , M w , m b ) 



-- 1 r~ b B^{ - s; M w , M w ) - r tw 



Bq{ - s;m b ,m b ) - 2 



B^(-m 2 ;m b ,M w ) - 1 
- (2 + r bw ) B% w ( - m 2 ; m b , M, 



(1 - r tw ) (2 + r tw ) ~ r bw (1 - 2r tw + r bw ) 



MlB 0p ( - m 2 ; m b , M w ) 



2r bw (l + r^j ~^- L ab (rn t , m b ,M w ) 
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— c 



w 



f 4 (1 - r tw ) + K bw 



4 + 



s 2 -c 2 



(4 + %) 



xMlC (-m 2 t ,-m 2 t ,-s;M w ,m b ,M, 



m t 2 ; M w , M w ) - - m 2 ; m 6 , M w ) 



B dw(~ m 2 ;m b ,M w ) + 1 



-2B^(-m 2 ;m b ,M w ) 



~U+ I2r tw - ^-^r tw (7 + r^) 

2 V c i 



1 



s 2 -c 2 
4+^4-^(1-tw) 



M 2 

(m t ,m b , M, 

^3r 



^3(Ti - 5 t Ajr Im ^ImS^( - m 2 ; M w , m b ) , 



(2.79) 



TI' 



4M^C ( - m 2 , -m 2 , -s; m 6 , M ff) m b ) 



M 2 1 

--Bjf( - s; m b , m b ) + B$( - m 2 ; m b , M w ) + 1 - 2 (3 + -^ L ab (m t , m b , M w ) J 



-- (2 + r bw ) B% w ( - m 2 ; m b , M H 
1 

~4 



(1 - r tw ) (2 + r t6 ) + r 6w r t6 M^ J B 0p ( - m 2 ; m 6 , Af w ) 

{2 2 
(rwM^ - m 2 , -m 2 , -s; M w , m b , M w ) 
c w ^ 



<( - s; M w , M w ) - <( - m 2 ; m 6 , M w ) + 1 

M 2 



s 2 -c 2 



(3 + r tw + r bw ) 



—L na (m t ,m b ,M w 



*3r 



+ ^ A ?iJ lmB o( - m h M w, m b ), 



(2.80) 



1 J Vb + «fe 
Asrl 2/ t (3) 



8 + ( r tt) 2 - 6r w + 



--(10 + r tw -3r bw ) 



Rw 



M^C ( - -m t 2 , s; m b , M w , m b ) 



-,(2 + rS) 



S ( - s; m b , ra 6 ) - B (-m t ; m b , M M 
- m 2 ;m b , M w ) + 1 



-3 



(1 + r tw ) (2 + r tw ) - r bw (1 + r 6w ) 



M 2 

-r^L ab (m t ,m b ,M w ) 

^3r 
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+r bw ^M^C ( - m%, -m 2 , -s; m b , M w , m b ) 
1 



+ - 



B Q ( — s; m b , m b ) - B ( — m t ; m b , M w ) - 1 



M 2 

-3 (l + r- b ) -r^L ab (m t , m b , M M 



-c,. 



4r bw - 



s 2 -c 2 



1 



,1 — r tw) (2 — r tw ) — r bw ( 5 — r tw — Ar bw — 



xMlC ( - m 2 , -m 2 , -s; M w , m b , M w ) 



1 

+ 2 



4 - '^J^ (4 - r m - 5r bw ) 



('- 



<( - s; M w , M w ) - B%{ - m 2 ; m b , M w ) 



if s 2 -c 2 

+ 5 4 "^ r 



B dw (-m 2 ;m b ,M W/ 



-6 1-r 



M 2 



A 



•"no 



3r 



Here we introduce more symbols, which were not given in Eq. ([2.10|) : 

r i& = r iw ± rw , A 3r = Am 2 t - s . 

Furthermore, we used one more 'subtracted' function: 

1 



(2.81: 



(2.82) 



5 L( -m 2 ;m b ,M w 



~ r bw 



1 - rw) [S F ( - m 2 - m b , M w ) + L M (M 2 .) - 1 
L M (m 2 ) - L,(M 2 W 



(2.83) 



and the three auxiliary functions: 

L ab (M 1} M 2 ,M 3 ) = (Ml + Ml -M 2 )c (-ml-m 2 } -s;M 2 ,M 3 ,M 2 ) 



-B* ( - s; M 2 , M 2 ) + fl* ( ~ m t 2 ; M 2 , M- 



3), 



L na (M u M 2 ,M 3 ) = (M 2 -M 2 -M 2 )C (-m 2 ,-m 2 ,- S] M 3 ,M 2 ,M 3 

+< ( - s; M 3 , M 3 ) - < ( - m 2 ; M 3 , M 2 ), 



L ffi (M!,M 2 ,M 3 ) 



- (M| + M 3 2 ) - 2M 2 



C ( - m 2 , -m 2 , -s- M 2 , M 1 , M 3 ) 



(2.84) 

(2.85) 
(2.86) 



+<( - s; M 3 , M 2 ) - ~B£( - ml M h M 2 ) - -B%{ - m 2 ; M 3 , M x ). 

Four scalar form factors, J-q[d , as follows from calculations, are both gauge-invariant and 
finite, thus enlarging the number of gauge-invariant building blocks to 15. On the contrary, 
two form factors, J-2^ w , are neither gauge-invariant nor finite. Gauge dependence on £, as 
well as UV poles, of L form factors cancel in the sum with the WW box and the self-energy 
contributions. 
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2.5 Library of scalar form factors for electron vertex 



Besides Btt clusters, we need also Bee clusters, which can, in principle, be taken from || or 
derived from the Btt case in the m t — > limit. Here we simply list the results: 



with 



-—Q e v e a e T z -*{s)+T w —{s) 



5 2 F Z " 



4cL 



2c^ 



A e ^ z - (s) + (s) , 



4c 2 

w 



In Eq. ( |2.87 ) we use three more auxiliary functions: 

\2 



J^ Z e EE J^ 2 



o (l + ^. 



(2.87) 
(2.88) 



-MjC (0, 0, -s; 0, M z , 0) - 3 <( - s; 0, 0) + L^(M 2 Z ) 



B%(-s;0,0) + L tM (M 2 z )-l 



[2.89) 



-2R V 



M 2 C (0, 0, -s; M w , 0, M w ) + B*{ - s; M w , M w ) + L^(M 2 



-4M^C (0, 0, -s; M w , 0, M w ) - S F ( - s; M w , M w ) - 3L M (M£) + - , (2.90) 



' /,V) -M^C o (0, 0, -s- 0, M w , 0) - \ [<( - s; 0, 0) + L,(M 2 W ) 



R v 

i-R 

4 Uw 



B^(-s;0,0) + L,(M 2 w )-l 



[2.91) 



2.6 The WW box 

There is only one, crossed, WW diagram contributing to our process, see Fig. | 
Here we give the contribution of this diagram to the scalar form factor LL: 



(2t) 4 i 9 1 

167T 2 S 



7^7+ ® 7/,7+^l 



where 



(w — rrclj -Do(0, 0, — m 2 , — m 2 , —s, —u; M w , 0, M w , m b ) 



+C ( - m 2 , -m 2 , -s; M w , 0, M w ) + C (0, 0, -s; M w ., 0, M w ) 



(2.92) 



(2.93) 
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b, s, d 



e W t 

Figure 9: Crossed WW box. 



with s, t, and u being the usual Mandelstamm variables satisfying 



s + t + u = 2m\ . 



2.7 The ZZ box 



(2.94) 



There are four ZZ diagrams, which form a gauge-invariant and UV finite subset. Its contri- 
bution is originally presented in terms of six structures (L, R) <g> (L, R, D) (i.e. here we used 
initially the L, R, D basis): 



B 1 



d+c 



1 «? 4 1 



32 c^s 



[7^7+ ® 7^7+] Fll 0, * , u) + [ 7m 7+ ® 7„7-] ^ («, * , «) 
+ [7^7- ® 7 M 7+] Kl («, * , «) + [7m7- ® 7^7-] ^ (s, *, «) (2.95) 
+ [7m7+ ® (-imtlD^)} TH (s, t, u) + [ 7m7 _ <g> (-im* (s, t, u) 




+ 



fWVWWW> 



► — *\aaa/vww\* — ► 

e " z * 



e z t 

Figure 10: Crossed and direct ZZ boxes. 

Moreover, we used three auxiliary functions J 7 , Ti, Q: 

Fll(s,t,u) = J r (a e ,a t ,5 t ,s,t)-J : '(a e ,5 t ,a t ,s,u), 



30 



F z L z L {s,t,u 



T (S e , S t , o t , s,t) - T (S e , a t , S t , s, u) , 
Ti (a e , a t , S t , s,t)-Q (a e , S t , a t , s, u) , 
H (S e , S t , o t , s,t)-Q (S e , a t , S t , s, u) , 
Q (a e , a t , S t , s,t) -H (a e , S t , a t , s, u) , 
Q {S e , S t , o t , s,t) -H (S e , o t , S t , s, u) . 



(2.96) 



Separating out Z fermion coupling constants and some common factors, we introduce more 
auxiliary functions. For T (<r e , o t , S t , s, t) defined as 



Mr 



F{<j e ,<j t ,St,s,t) = — o e a t T x (s,t) + a e a t StT 2 (s,t) 



(2.97) 



there are two 
Fi(s,t) = 



t + M 2 + t 2 _ 



-s z + 2M 2 ) t 2 _ + AtM 



t- 



A 



4r 



xD Q {0, 0, -m 2 , -m 2 , -s, -t; M z , 0, M z , m t ) 



" 2 M 2 t (t-t + 2M 2 t + ) + sm\ , 2t + (M 2 z -2m 2 t ) 



A 



+ 



4r 



*3r 



xC ( - m 2 , -to 2 , s; M z , m t , M z ) 



t 1 + 



Ml 



+ 



t 3 + 2tt_M 2 - m 2 (Mt- + m\) 



C (0,-m 2 ,t;0,M z ,m t ) 



9 (t - 2t_ - 2M 2 ) t 2 + smf 
Ml - to? + z ' - 



z ""t 

s + 2t_ 



i4 r 



C (0,0,s;M z ,0,M z 



1 + 



^3r 



5 ( - s; M z , M z ) - B ( - m 2 ; M z , m t ) 



2t 



B o (t;m t ,0) -B (-m 2 ;M z ,m t ) 



(2.98) 



and 



T 2 {s,t) = -{t 2 ^ + 2M 2 z t)D {0,0,-m 2 ,-m 2 ,-s,-t;M z ,0,M z ,m t ) 
-t+C ( - m 2 t , -m\, s; M z , m t , M z ) 
+2tC (0, -m 2 , t; 0, M z , m t ) - t_C {0, 0, s; M z , 0, M z ). 

For H written out as 



(2.99) 



H(a e ,a t ,St,s,t) 



x 4r 



aiapix (s, t) + a 2 e a t S t H 2 (s, t) + a 2 e S 2 H 3 (s, t) , (2.100) 



we need three auxiliary subfunctions: 
~st- 



Hi (s,t) 



- - (t_ + MJ) 2 - 2M|t - 
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xD o (0, 0, -m 2 t , -m 2 , -s, -t; M z , 0, M z , m t ) 



t 



tl + A 3r m 



+ 



2A 



\C (-ml -m\ ,s;M z ,m t ,M z ) 



lr 



+ (2t + t- + t 



2M|t ^ + V o(0)0;S;Mz)0)Mz) 



Mr 



*- + 



(sm 2 — 



2A 



lr 



C (0,-mlt;0,M z ,m t ) 



(2.101) 



771 / \ 

-2j±- 1 £„(*; m t , 0) - B ( - ml M z , m t ) J - B (*; m t , 0) + B (s; M z , M z ), 



n 2 (s,t) 



f_ (s + t_) + 2m 2 M 2 



D (0, 0, -m?, -m t 2 , -s, -t; M z , 0, M z , m t ) 
- (s + t_ - 2m 2 ) C ( - m t 2 , -m 2 , s; M z , m t , M z ) (2.102) 
-2m 2 C ( - ml -ml s; 0, M z , m t ) + (s + C o (0, 0, s; M z , 0, M z ), 



and 



n 3 (s,t) 



—s 



t-D o (0, 0, -ml -ml -s, -t; M z , 0, M z , m t ) 



+C ( - ml -ml s; M z , m t , M z ) + C (0, -m 2 , t; M z , 0, M z ) 



Finally, (/ also, denned as follows: 

Q (a e ,a t ,8 t ,s,t) 



needs only two additional functions: 



\4 r 



a e a t G 1 (s,t) + a e a t 5 t Q 2 (s, t) 



(2.103) 



(2.104) 



Qi (s,t) 



,2 2 



- tA + tM 2 (2m 2 - M 2 ) 



xL> (0, 0, -ml -ml -s, -t; M z , 0, M z , m t ) 



+ 



t + (t + t_) + 2M> 2 f 



C ( - ml -ml s; M z , m t , M z ) 

C (0,-ml-t;0,M z ,m t ) 



s z ( t + t 2 

t-t + —[t- -±— 

2 V A 4r 



q,(0,0,s;M z ,0,M z ) 



+2mf 1 + 



t- 



B (t; m t , 0) - B ( - m 2 ; M z , m t ) 



B o (s;M z ,M z )-B o (t;m t ,0) 



(2.105) 
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and 



02 (M) 



rn 



(t 2 _ + 2M 2 z t)D (0, 0, -m 2 , -m 2 , -s, -t; M z , 0, M z , m t ) 
+t+C ( - m t 2 , -m 2 , s; M z , m £ , M z ) - 2tC (0, -m 2 , £; 0, M z , m*) 



+t_C (0, 0, s; M z , 0, M z 
In this section we used the notation: 



Mr 



(2.106) 



(2.107) 



together with of Eq. (|2.82| ), this denotes remnants of Gram determinants that remained 
after cancellation of factors s and 4, leading to a simplification of the expressions. 

2.7.1 Transition to the L, Q,D basis 

Since the ZZ box contribution is given in the L, R, D basis, while all the rest is in the L, Q, D 
basis, we should transfer one of them to a chosen basis. At this phase of the calculations there 
is not much difference which basis is chosen. For definiteness we choose the L, Q, D basis and 
transfer the ZZ box contribution to it. The transition formulae are simple: 



■pzz 

LL 


{*, 


t, 


u) 




t, 


u 


) - 


F z L z (s,t,u)-F z R z L (s,t,u 


'pzz 

QL 


(«, 


t, 


u) 


= 2 


3~ RL ( S 5 ^5 U ) 3~ RR 




t, 


uj 


1 


<pzz 

LQ 




t, 


u) 


= 2 


3~LR ( S 5 ^1 U ) ~ 3~ RR 




t, 


u) 




'pzz 

•'qq 


(s, 


t, 


u) 










■ttzz 

LD 


(*, 


t, 


u) 


3~LD ( S 5 ^1 U ) J~ RD i S - 


t, 


u 


), 


rzz 

QD 


(«1 


t, 


u) 











(2.108 
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3 Scalar form factors for electroweak amplitude 

Having all the building blocks, it is time to construct complete electroweak scalar form factors. 



3.1 Vertices scalar form factors 

We begin with two vertex contributions: 




Figure 11: Electron (a) and final fermion (b) vertices in ee — > (7, Z) — > //. 

In the same way as described in ||, we reduce two vertex contributions to our six form 
factors: 



F LL (s) 



F QL (s) 



F QQ 00 



Fr(s) + F z L tt (s)-4c 2 w A(M, 



F z Q ee (s) + Ff (s) - 2c 2 w A(M w ) + k 



F? tt 



2A(M M 



F z L ee (s) - 2c 2 w A(M w ) + Ff (s) + k [F2 ee (s) - 2A(M W )} 

k 



is) + Ff (s) 



F£° {s) + Ff (s) 



where 



F LD (s) 



F QD (s) 



jiZtt 
D 



Ff (s) + kF2 tt (s) 



(3.109) 



k = c 2 w (R z - 1) . 
With the term containing A(M W ), 



A(M M 



1 M 2 
4-ln^f 

e \i l 



(3.110) 



(3.111) 



we explicitly show the contribution of the so-called special vertices |^TJ . Note that they ac- 
company every L form factor. The poles 1/e originating from special vertices will be canceled 
in the sum of all contributions, including self-energies and boxes. 
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3.2 Bosonic self-energies and bosonic counterterms 

The contributions to form factors from bosonic self-energy diagrams and counterterms, orig- 
inating from bosonic self-energy diagrams, come from four classes of diagrams; their sum is 
depicted by a black circle in Fig. |TJ. 





+ 





+ 




Figure 12: Bosonic self- energies and bosonic counterterms for ee — > (Z,j) — > ff. 



The contribution of these diagrams to the four scalar form factors is derived straightfor- 
wardly 0, §: 



Ft (*) 



rpct 
QL(LQ) 



r>rf,bos 



QQ 



net. fer 



c 2 -s 2 



V z (s) - s 2 w U J7 (0) + ., 



Ap + Ap bos ) 



v z (s)-(n zi (s) + n h z °s 



.s)i -4n 77 (o)-(A P + Ap bos 



v h ™ (s) - 2 (n b ° s ( s ) + n bo ; ( S )) + k[u^ s ( s ) - n b 7 s (o) 
-4.n b 7 s (o) - -1 (Ap bos + Ap 



—bos 



(s) - 2U%(s) - sitf- (0) - — Ap fcr . 



(3.112) 
(3.113) 

(3.114) 
(3.115) 
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We note that the term (s) — ife (s)} is conventionally extracted from F^ iei (s). This 



contribution is shifted to A™ A , Eq. ( |1.8| ). 

In Eqs. ( |3.112| )-( p.ll5| ) Ap bos and fi^ 08 (s) stand for shifts of bosonic self-energies. They 
have the same origin as special vertices and they are equal to: 



bos 



As 2 w A(M w ), (3.116) 
-2R W A(M W ), (3.117) 

see Eqs. (6.137) and (6.139) of 0. These poles also cancel in the sum of all contributions. 



n bos (s) 

z~ f \ / 



3.3 Complete scalar form factors of the one-loop amplitude 

Adding all contributions together, we observe the cancellation of all poles. The ultraviolet- 
finite results for six scalar form factors are: 



F LL (s,t,u) 



Fql (s, t, u) 
F LQ (s, t, u) 



TT («) + Tf (s) + Tt (s) + k™ K L W (s, u) + k zz TH (s, t, u) , 

TT ( s ) + J** ( s ) + k jr T e {s) + {s) + k zz T zz (s> tf u) ) 



F QQ ( S > U ) 



Tr(s) + T z n tt (s)- — 



u ( 



"yee 



F LD (s, t, u) 



F QD (s, t, u) 



Tf{s) + kTf{s) + k zz F z Q z D {s,t,u) 



(3.118) 



where 



16k. 



(3.119) 



k 1 



2cl 



(3.120) 



In Eq. ( |3.118|) , the quantities T^ B (s) denote finite parts of the counterterm contributions, see 
Eqs. ( PT2D - (|3TTT5|) . 

The formulae of Sections 2 and 3 put together present the one- loop core of the eef fLib 
code. 
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4 Improved Born Approximation cross-section 

In this section we give the improved Born approximation (IBA) differential in the scattering 
angle cross- sect ion. It is derived by simple squaring the (7 + Z) exchange IBA amplitude, 
Eqs. (|1.8|) - (|TTniD , and accounting for proper normalization factors. We simply give the result: 



da 1 



dcos'd 



Aiv c « A + <! A + v 



IBA\ 
ZZ ) 1 



(4.121) 



a 



where — y 1 — Am 2 / s and 

Q 2 (s 2 + 2st + 2t 2 _) \a (s) 
2Q t Rej X (2 L 

+ (s 2 + 2st + 2t 2 _ 



a 



IBA 

77 

JBA 



7 Z 



s + t-) 2 + sm 2 



F LL (s,t,u) 



(s, t, u) + F L q (s, t, u) + Fqq (s, t, u) 



Am 2 (st + t: 



F LD (s,t,u) + Fq D (s,t,u) 



a is. 



\ X \ 2 Re{8 (s + t 
+2 



s + t) 2 + t 2 _ 



F LL (s,t,u] 



Fql (S,t,u) 



+ F LL (s,t,u)F* (s,t,u) 



+4 



[s + t-) + sm 2 



2Fr. r .(s,t,u)F*(8,t,u) 



LQ ") 

+F LL (s, t, u) F* Q (s, t, u) + F QL (s, t, u) F*q (s, t, u] 

2 



+ 



2 st + t 2 



F LQ (s,t,u) 



+ 



Fqq (s, t, U / 



+2 F QL (s, t, u) + F L q (s, t, u) F* (s, t, u 



-8m 2 (st + 1 2 



2F LD (s, t, u) + Fq D (s, t, u) F* LL (s, t, u) 



+ [F LD (s,t,u) + Fq D (s,t,u) IP* (s,t,u 



+ ( 2F LD (s, t, u) + F QD (s, t, u) ) Ft n (s, t, u 



LQ "1 



+ F LD (a, t, u) + Fq D (s, t, u) ) F* (s, t, u 



-2m 2 (st + t 2 _ ) A 3r 



F LD (s, t, u) 



+2F LD (s, t, u) F* (s, t, u) + F QD (s, t, u) 



(4.122) 
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5 Numerical results and discussion 



All the formulae derived in this article are realized in a FORTRAN code with a tentative name 
eef f Lib. All the numbers are produced with December 2000 version of the code [pl|. In this 
section we present several examples of numerical results. 

We will show several examples of comparison with ZFITTER v6.30 [19]. In the present 
realization, eef f Lib does not calculate M w from \x decay and does not precompute either 
Sirlin's parameter Ar or total Z width, which enters the Z boson propagator. For this reason, 
the three parameters: M w , Ar , V z were being taken from ZFITTER and used as INPUT for 
eef f Lib. Moreover, present eef f Lib is a purely one- loop code, while in ZFITTER it was 
not foreseen to access one-loop form factors with users flags. To accomplish the goals of 
comparison at the one-loop level, we had to modify a little the DIZET electroweak library. The 
most important change was an addition to the SUBROUTINE ROKANC: 
* 

* For eett 
* 

FLL=(XR0K(1)-1D0+DR )*R1/AL4PI 
FQL=FLL+ (XROK (2) -1D0) *R1/AL4PI 
FLQ=FLL+ (XROK (3) -1D0) *R1/AL4PI 
FQQ=FLL+ (XROK (4) -1D0) *R1/AL4PI 

with the aid of which we reconstruct four form factors from ZFITTER's effective couplings p 
and k's (F ld and F QD do not contribute in massless approximation). 



5.1 Flags of eeffLib 

Here we give a very brief description of flags (user options) of eeffLib. While creating the 
code, we followed the principle to preserve as much as possible the meaning of flags as described 
in the ZFITTER description ||. In the list below, a comment 'as in ZFD' means that the flag 
has exactly the same meaning as in ||. 



ALEM=3 ! as in ZFD 
ALE2=3 ! as in ZFD 
VP0L=0 

Note that the flag is extended to VP0L=0 to allow calculations 'without running of a\ 



=0 \alpha(0) ; =1,=2 as in ZFD; =3 is reserved for later use 



QCDC=0 
IT0P=1 
GAMS=1 
WEAK=1 
IM0MS=1 



as in ZFD 

as in DIZET (internal flag) 
as in ZFD 

as in ZFD (use WEAK=2 in v6.30 to throw away some H0-terms) 
=0 \alpha-scheme ; =1 GFermi-scheme 



New meaning of an old flag: switches between two renormalization schemes; 
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BOXD=0 ! =0 without any boxes 

=1 with 77 box 

=2 with Z^i box 

=3 with 77 and Z7 boxes; 1, 2, 3 are used together with WEAK=0 

=4 with WW box 

=5 with WW and ZZ boxes; 4, 5 are used together with WEAK=1 

GAMZTR=1! =0 GAMZ=0; =1 GAMZ.NE.O 

Treatment of Y z . The option is implemented for the sake of comparison with FeynArts. 

EWFFTR=0! =0 EWFFs ; =1 RHO-KAPPAS 

Treatment of EW form factors; switches between form factors and effective ZFITTER 
couplings p and «'s. The option is implemented for comparison with ZFITTER. 

FERMTR=1! =1 a 'standard' set of fermions masses; =2,3 'modified' 

Treatment of fermionic masses; switches between three different sets of 'effective quark 
masses'. 



Table 1: EWFF for the process e+e~ — >■ uu. eef f Lib-ZFITTER comparison. 



Without EW boxes 


Quantity 


E 

- LJ cm 




100 GeV 


200 GeV 


300 GeV 


F LL 


M w /10 
M w 
10M W 


13.47777 ■ - zl.84781 
13.47777 - il.84781 
13.47777 - il.84781 


16.22034 -ilO.49408 
16.22034 -ilO.49408 
16.22034 -ilO.49408 


23.75241 -ill.27464 
23.75241 -ill.27464 
23.75241 - ill.27464 


ZFITTER 


13.47771 - il.84786 


16.22031 -ilO.49405 


23.75237 -ill.27464 


Fql 


M w /10 
M w 
WM W 


29.34725 + i3.67334 
29.34725 + i3.67334 
29.34725 + i3.67334 


30.33892 + i3.34535 
30.33891 + i3.34535 
30.33891 + i3.34535 


31.64554 + %2. 75260 
31.64554 + %2. 75260 
31.64554 + %2. 75260 


ZFr 


rTER 


29.34720 + i3.67330 


30.33889 + i3.34535 


31.64552 + i2. 75259 


Flq 


M w /10 
M w 
10M W 


29.13302 + i3.26972 
29.13302 + i3.26972 
29.13302 + i3.26972 


30.03854 + il.54158 
30.03854 + il.54158 
30.03854 + il.54158 


31.68636 - i0.22635 
31.68636 - i0.22635 
31.68636 - i0.22635 


ZFITTER 


29.13304 + i3.26973 


30.03855 + il.54163 


31.68635 - i0.22634 


Fqq 


M w /10 
M w 
WM W 


44.90390 + i8.85688 

44.90389 + i8.85688 

44.90390 + i8.85688 


43.80287 +il0.02412 

43.80285 +il0.02412 

43.80286 +il0.02412 


44.21224 + il0.83899 

44.21222 + il0.83899 

44.21223 + il0.83899 


ZFr 


rTER 


44.90392 +i8.85688 


43.80285 + ilO.02411 


44.21224 + il0.83894 


WW is added 


F LL 


M w /10 
M w 
10M W 


12.94471 - il.84781 
12.94471 - il.84781 
12.94471 - il.84781 


9.34003 - i9.42493 
9.34003 - i9.42493 
9.34003 - i9.42493 


9.03774 -ill.56004 
9.03774 -ill.56004 
9.03774 -ill.56004 


ZFITTER 


12.94468 - il.84786 


9.34065 -- i9.42467 


9.03903 -ill.55958 
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Table 2: EWFF for the process e + e — > uu. eef f Lib-ZFITTER comparison. 



With ZZ boxes 


Quantity 


F C m 




100 GeV 


200 GeV 


300 GeV 


F LL 


M w /10 
M w 
10M W 


12.89587 - z'1.84781 

^ r > RQ^Rf, — i'1 84781 

12 89587 - z'l 84781 


8.24674 -z'10.64677 
8 9zLfi7^ 7'1DfizLfi77 

8 24673 - z'10 64677 


8.98241 -il2.88512 

8 Q89/L1 i'1 9 88^1 9 

8 98241 - H2 88512 


ZFITTER 


12.89583 - z'1.84786 


8.24736 -ilO.64651 


8.98370 -il2.88466 


Fql 


M w /10 
M w 
WM W 


29.30451 + z3.67334 
29.30451 + i3.67334 


29.38219 + z'2. 27613 

9Q^8918 J- i9 97fi1? 

29.38219 + z'2. 27613 


31.59712 + z'1.59304 
31.59712 + z'1.59304 


ZFITTER 


29.30445 + z'3.67330 


29.38216 + z'2. 27613 


31.59710 + z'1.59304 


Flq 


M w /10 
M w 
10M W 


29.10829 + i3.26972 
29.10829 + z3.26972 
29.10829 + i3.26972 


29.48510 + z0.92306 
29.48509 + i0.92306 
29.48509 + z0.92306 


31.65836 - i0.89713 
31.65835 - i0.89713 
31.65835 - i0.89713 


ZFITTER 


29.10832 + i3.26973 


29.48512 + i0.92312 


31.65835 - z'0.89711 


F QQ 


M w /10 
M w 
10M W 


44.88226 + i8.85688 
44.88226 + i8.85688 
44.88226 + i8.85688 


43.31856 + z9.48287 

43.31854 + i9.48287 

43.31855 + z'9.48287 


44.18773 +Z10.25200 

44.18771 +Z10.25200 

44.18772 +z!0.25200 


ZFITTER 


44.88228 + i8.85688 


43.31854 + i9.48286 


44.18773 +Z10.25196 



5.2 eef f Lib— ZFITTER comparison of scalar form factors 

First of all we discuss the results of a computation of the four scalar form factors, 

F LL (s,t), F QL (s,t), F LQ (s,t), F QQ (s,t), (5.123) 

for three variants: 

1) without EW boxes, i.e. without gauge-invariant contribution of ZZ boxes, and without 
f = 1 part of the WW box, Eq. fl2J3|); 

2) without ZZ boxes; 

3) with full content of EWRC. 

In this comparison we use flags as in subsection and, moreover, 

M w = 80.4514958 GeV, 
Ar = 0.0284190602, 

T z = 2.499776 GeV. (5.124) 

In Table || we show an example of comparison of four form factors F LL QL LQ QQ (s, t) between 
the eeffLib, where we set m t = 0.2 GeV and ZFITTER (the latter is able to deliver only 
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massless results). The form factors are shown as complex numbers for the three c.m.s. energies 
(for t = mf — s/2) and for the three values of scale \i = M w /10, M w , 10M W . The table 
demonstrates scale independence and very good agreement with ZFITTER results (6 or 7 digits). 
One should stress that total agreement with ZFITTER is not expected because in the eef f Lib 
code we use massive expressions to compute the nearly massless case. Certain numerical 
cancellations leading to losing some numerical precision are expected. We should conclude 
that the agreement is very good and uniquely demonstrates that our formulae have the correct 
m t — > limit. 

In Table |2| we show a similar comparison with ZFITTER when ZZ boxes are added. As 
seen, the agreement has not deteriorated. 



5.3 eef f Lib— ZFITTER comparison of IBA cross-section 

As the next step of the comparison of eef f Lib with calculations from the literature, we present 
a comparison of the IBA cross-section. 

In Table ^] we show the differential cross-section Eq. ( |4.121| ) in pb for three values of cos d = 
—0.9, 0, +0.9, with IPS of Eq. ( |5.124j ) and without running e.m. coupling, i.e. a (s) — > a. 



Table 3: IBA, First row - ZFITTER (uu channel); second row - eef f Lib (m t = 0.1 GeV); third 
row - eef f Lib (m t = 173.8 GeV). 





lOOGeV 


200GeV 


300GeV 


400GeV 


700GeV 


lOOOGey 


cos ?9 = -0.9 


47.664652 


0.291823 


0.169510 








47.661401 


0.291827 


0.169515 


0.103284 


0.035318 


0.017203 








0.162579 


0.043974 


0.018850 


cos $ = 


59.768387 


1.718830 


0.695061 








59.770715 


1.718870 


0.695072 


0.376868 


0.117276 


0.055870 








0.264874 


0.112923 


0.054211 


costf = 0.9 


168.981978 


5.954048 


2.292260 








168.991272 


5.954166 


2.292289 


1.222343 


0.372903 


0.176030 








0.438952 


0.293415 


0.154784 



Next, we present the same comparison as in Table |3|, but now with running e.m. coupling. 
Since the flags setting VP0L=1, which is relevant to this case, affects ZFITTER numbers, we 
now use, instead of Eq. Q5.124Q , the new INPUT set: 

M w = 80.4467671 GeV, 
Ar = 0.0284495385, 

T z = 2.499538 GeV. (5.125) 

The numbers, collected in Table |], exhibit good level of agreement. 

Finally, in Table we give a comparison of the cross-section integrated within the angular 
interval | cos$| < 0.999. (Flags setting is the same as for Table 4.) 
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Table 4: IBA, First row - ZFITTER (uu channel); second row - eef f Lib (m t = 0.1 GeV); third 
row - eef f Lib (m t = 173.8 GeV). 





lOOGeV 


200GeV 


300GW 


mGeV 


imGeV 


lOOOGe^ 


cos i9 = -0.9 


A K AC\A7 AO 

40.4U4 ( 4z 
45.404593 


U.ooDyoo 
0.386966 


0.225923 


0.138065 
0.195069 


0.048621 
0.057892 


0.024155 
0.025877 


cos-# = 


60.382423 
60.382553 


1.882835 
1.882835 


0.771939 
0.771938 


0.421409 
0.303984 


0.133474 
0.130208 


0.064244 
0.062853 


cos$ = 0.9 


173.467517 
173.467515 


6.450000 
6.449995 


2.510881 
2.510877 


1.346616 
0.493006 


0.417292 
0.330480 


0.198839 
0.175598 



Table 5: eef f Lib-ZFITTER comparison of the total cross-section. Cross-sections are given 
in picobarns: the first row - aH t , i.e. eef f Lib (mt = 0.1 GeV); the second row - cr^f, i.e. 
ZFITTER (mm channel); the last entry shows the deviation (al[ t — er^f J /cr^f in per mill. 



100 GeV 


200 GeV 


300 GeV 


Otot 


0"fb 


Otot 




Ctot 


0"fb 


160.8981 
160.8980 


70.98419 
70.98406 


5.021797 
5.021808 


3.360836 
3.360848 


2.031750 
2.031754 


1.269552 
1.269556 


+0.001 


+0.002 


-0.002 


-0.004 


-0.002 


-0.003 



A typical deviation between eeffLib and ZFITTER is of the order ~ 10 -6 , i.e. of the 
order of the required precision of the numerical integration over cos$. Examples of numbers 
obtained with eeffLib, which were shown in this section, demonstrate that ZFITTER numbers 
are recovered for light m t . 

We conclude this subsection with a comment about technical precision of our calculations 



(modulo bugs, of course). We do not use looptools package [23]. For all PV functions, but 



one, namely Dq function, we use our own coding where we can control precision internally 
and, typically, we can guarantee 9-10 digits precision. For D function we use, instead, REAL*8 
T0PAZ0 coding |22j and the only accessible for us way to control the precision is to compare 
results with those computed with REAL* 16 T0PAZ0 coding. This was done for a typical Do 
function entering WW box contribution. We got an agreement within 9 digits between these 
two versions for all = 400, 700, 1000 GeV and cos# = 0.9, 0, -0.9. 



5.4 About a comparison with the other codes 

As is well known, the one-loop differential cross-section of e + e~ — > ti may be generated with 
the aid of the FeynArts system |23| . FeynArts-generated versions of the code with and without 
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QED contributions are available and an attempt to compare results was undertaken. We 
compared da/dcos8 without QED contributions at yfs = 400, 700, 1000 GeV and three values 
of cos 6 = 0.9, ,0, —0.9. An agreement of numbers for the Born cross-section within 6 digits 
was found, while for the one-loop corrected cross-section we managed to reach an agreement 
within 1 — 3% only. 

There is another FORTRAN code for ti production. It was originally written for the MSSM 
25|j , but it has also been tailored for the SM. So far we managed to completely agree with this 



code only for the Born cross-section; while we turned to the one-loop one, we realized that no 
separation between QED and EW corrections is implemented into this code. On the contrary, 
in the eef fLib version used to produce numbers for this paper, we coded only the EW part 
of the cross-section. A present, the QED part is also available in our code |T7||. Moreover, this 



code produces cross-section integrated over the angle, while it would be more informative to 
compare the differential quantity: 

rf /— da (1) {s,t)/dt , uo „, 

^,cos^)= dai0){s)/dt "I- (5-126) 

For the time being we limit ourselves by presenting eef fLib results for <5(i/i, cos-$) of 
Eq. ( |5.126| ), which are shown in Fig. [T3|. (Flags setting is the same as for Table 3.) 

Furthermore, in Fig. 3 of paper |fL6| , an interesting result is presented. We tried to 
reproduce it with the aid of eef fLib. The results are shown in Fig. [14|. As might be seen from 
a comparison of two figures, there is nearly ideal agreement for in the interval [500-3000] 
GeV, while above 3000 GeV the eef fLib curve goes a bit higher than the curve shown in Hl6fl . 
Note, that both curves show a very similar M H dependence. It is difficult to expect more from 
such a pilot comparison, because even input parameters and various options were not tuned. 
Meantime, a Bielefeld-Zeuthen team P3| started alternative calculations using the DIANA 



system |[27|| . A comparison of results was undertaken. It showed good agreement of numbers. 



Recently, we were provided with the numbers computed with the FeynArts system |28| for 
da/d cos 9 without QED contributions; they showed better agreement than we managed to 
reach ourselves. 

The results of latest comparisons will be presented in more detail elsewhere. 
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Figure 13: Relative EWRC 5 [Eq. (|5.126|) 1 to the e + e — ► tt differential cross-section. Numbers 
near the curves show yfs in GeV. 
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Figure 14: Relative EWRC to e+e" -> for M ff = 100 GeV (solid line) and M H = 1000 GeV 
(dashed line). 
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